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Abstract. In the SK analysis of the neutrino events for Fully Contained Events and Partially Contained 
Events on their zenith angle distribution, it is assumed that the zenith angle of the incident neutrino is 
the same as that of the detected charged lepton. In the present paper, we examine the validity of the SK 
assumption on the direction of the incident neutrinos. Concretely speaking, we analyze muon-like events due 
to QEL. For the purpose, we develop Time Sequential Monte Carlo Simulation to extract the conclusion on 
the validity of the SK assumption. In our Time Sequential Simulation, we simulate every physical process 
concerned as exactly as possible without any approximation. From the comparison between the zenith 
angle distributon of the emitted muons under the SK assumption on the direction and the corresponding 
one obtained under our Time Sequential Simulation, it is concluded that the measurement of the direction 
of the incident neutrino for the neutrino events occurring inside the detector in the SK analysis turns out 
to be unreliable, which holds irrespective of the existence and/or non-existence of the neutrino oscillation. 

PACS. Superkamiokande, QEL, Fully Contained Event, Time Sequential Monte Carlo Simulation 



1 Introduction 

Superkamiokande have been analyzing Fully Contained 
Events and Partially Contained Events which are gen- 
erated inside the detector, and Upward Through Going 
Events and Stopping Events which are generated outside 
the detector, for the studies on the neutrino oscillation in 
atmospheric neutrinos. The report of oscillations between 
muon and tau neutrinos for atmospheric neutrinos de- 
tected with SuperKamiokande (SK, hereafter) is claimed 
to be robustly established for the following reasons: 

(1) The discrimination between electrons and muons in the 
SK energy range, say, several hundred MeV to several 
GeV, has been proved to be almost perfect, as demon- 
strated by calibration using accelerator beams pQ 0. 



1 The SK discrimination procedure between muon and elec- 
tron is constructed on the average value theory. In our opinion, 
discrimination procedure should be examined, taking into ac- 
count the stochastic characters of the physical processes in the 
neutrino events concerned. If we take this effect into account, 
then, for example, we give uncertainties of 3° to 14° in the in- 



(2) The analysis for the electron-like events and the muon- 
like events which give the single-ring structure in Fully 
Contained Events and Partially Contained Events with 
their zenith angle distribution, based on the well es- 
tablished discrimination procedure mentioned in (1), 
reveals a significant deficit of muon-like events but 
the expected level of electron-like events. It is, thus, 
concluded that muon neutrinos oscillate into tau neu- 
trinos which cannot be detected due to the small ge- 
ometry of SK. As the most new one, the SK collab- 
oration published their comprehensive paper [2]. The 
analysis of SK data presently yields sin 2 2# > 0.92 and 
1.5 x 1CT 3 cV 2 < Am 2 < 3.4 x KT 3 eV 2 at 90% confi- 
dence level. 

(3) The analysis of Upward Through Going Events and 
Stopping Events, in which the neutrino interactions 
occur outside the detector, leads to similar results to 

cident direction of the charged lepton and uncertainties of 2m 
to 7m in the vertex point of the events. See, our papers [2]. 
However, SK give 1.8° to 3.0° and 0.3m for the same physical 
quantities. See, accompanied two papers. 
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(2). The charged leptons which are produced in these 
categories are regarded as being exclusively muons, be- 
cause electrons have negligible probabilities to produce 
such events as they lose energy very rapidly in the sor- 
rounding rock. Thus for these events the discrimina- 
tion procedure described in (1) is not required, and, 
therefore, the analysis here is independent of the anal- 
ysis in (2). For these events, however, the SK group 
obtains the same parameters for neutrino oscillations 
as in (2) 0. 

(4) Now, SK assert that they have found the oscillatory 
signature in atomospheric neutrinos from L/E anal- 
ysis, which should be the ultimate evidence for the 
exsistence of the neutrino oscillation [3] Our critical 
examination on the L/E analysis by SK will be pub- 
lished elsewhere. 

As for item(3), we have clarified that SK hardly dis- 
criminate electron( neutrino) from muon( neutrino ) in 
the SK manner and, instead, propose more rigorous and 
suitable procedure with theoretical background for the dis- 
crimination between them in the preceeding two papers. 
Among the neutrino events both occurring inside and out- 
side the detector, the most robust evidence for the neu- 
trino oscillation, if exists, should have been obtained from 
the analysis of both electron-like events and muon-like 
events in Fully Contained Events. Because (i) all necces- 
sary informations for the physical interpretation are in- 
cluded in Fully Contained Events due to their character 
and (ii), 

furthermore, both electron-like events and muon-like events 
give the single structure image free from arbitrary inter- 
pretation with the proper electron/muon discrimination 
procedure. SK treat neutrino events whose energies cover 
from several hundred MeV to several GeV, if the neutrino 
events occur inside the detector. In this energy region, 
Quasi Elastic Scattering(QEL) [4] is dominant compared 
with other physical processes, such as one-pion produc- 
tion [5], coherent pion production [6] and deep inelastic 
scattering [7]. Events due to other processes, except QEL, 
are not free from ambiguities due to multi-ring structure 
of the images. 

Therefore, SK should have analyzed the muon- 
like events and the electron-like events with the 
single ring image in Fully Contained Events exclu- 
sively where QEL is dominant, without utilizing 
poorer quality events, if SK pursue to obtain the 
clear cut conclusion on the neutrino oscillation 0. 

2 It seems strange that the experimental data with different 
qualities give similarly precise results, because Fully Contained 
Events whose information is totally inside the detector are of 
higher experimental qualities compared with those of both Up- 
ward Through Going Events and Stopping Events. 

3 In their analysis, they really add Partially Contained 
Events to the experimental data as those with the same quality 
under the assumption that they belong to muon-like events in 
Fully Contained Events to raise the statistics higher. However, 
such the assumption lacks in theoretical background. Further- 
more, SK utilize to add multi-ring structure events which are 



Therefore, it is essential for us to examine single ring 
structure events among Fully Contained Events due to 
QEL which have the least ambiguities among the neutrino 
events concerned to obtain clear cut conclusion as for the 
neutrino oscillation. Here, the main concern of the present 
paper is devoted to the detailed analysis of the muon-like 
events from QEL, focusing on the direction of the incident 
neutrino. Situation around the corresponding electron-like 
event is the same as in the muon-like event. The exami- 
nation on the separation of Fully Contained Events from 
Partially Contained Events will be discussed in subsequent 
papers. 

Here, it should be emphasized that the direction of 
the incident neutrino is assumed to be the same as that 
of the emitted charged lepton, i.e., the (anti-)muon or 
(anti-) electron, in the SK analysis of both Fully Contained 
Events and Partially Contained Events 8,9J. The SK De- 
tector Simulation is to be constructed without any con- 
tradiction with the SK assumption on the direction. 

From the point of orthodoxical Monte Carlo Simula- 
tion, it seems to be unnatural for SK to impose such the 
assumption that the direction of the incident neutrino is 
the same as that of the emitted lepton ( hereafter, we call 
this assumption simply " the SK assumption on the di- 
rection" ) upon their Detector Simulation. Obviously, one 
need not any assumption on the relation between the di- 
rection of the incident neutrino and that of the emitted 
lepton in any sense, if we develop the Monte Carlo Simu- 
lation in a rigorous manner, which will be shown later in 
the present paper. 

In order to avoid any misunderstanding toward the SK 
assumption on the direction we reproduce this assumption 
from the original SK paper: 

" However, the direction of the neutrino must be estimated 
from the reconstructed direction of the products of the neu- 
trino interaction. In water Cherenkov detectors, the direc- 
tion of an observed lepton is assumed to be the direction of 
the neutrino. Fig \ll \ and Fia \12\ show the estimated corre- 
lation angle between neutrinos and leptons as a function of 
lepton momentum. At energies below 400 MeV/c, the lep- 
ton direction has little correlation with the neutrino direc- 
tion. The correlation angle becomes smaller with increas- 
ing lepton momentum. Therefore, the zenith angle depen- 
dence of the flux as a consequence of neutrino oscillation 
is largely washed out below 400 MeV /c lepton momentum. 
With increasing momentum, the effect can be seen more 
clearly. " 00. 

On the other hand, Ishitsuka states in his Ph.D thesis 
which is exclusively devoted into the L/E analysis of the 



caused by one-pion roduction, coherent pion production and 
deep inelastic scattering. However, the discrimination among 
the multi-ring structures is not so easy, which may lead the 
worse estimation of energies as well as directions of the events 
concerned. 

4 It could be understood from this statement that SK justify 
the validity of this assumption above 400 MeV/c. However, it is 
not correct, because SK put "to be proved " as the proposition. 
See, page 101 in their paper [8]. 
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atmospheric neutrino from Super Kamiokande as follows: 

" 8.4 Reconstruction of L v 

Flight length of neutrino is determined from the neu- 
trino incident zenith angle, although the energy and the 
flavor are also involved. First, the direction of neutrino is 
estimated for each sample by a different way. Then, the 
neutrino flight lenght is calclulated from the zenith angle 
of the reconstructed direction. 

8.4-1 Reconstruction of Neutrino Direction 
FC Single-ring Sample 

The direction of neutrino for FC single-ring sample is 
simply assumed to be the same as the reconstructed direc- 
tion of muon. Zenith angle of neutrino is reconstructed as 
follows: 



cos 6>r, ec = cos 0,. 



3.17) 



,where cos0£ ec and cos 0^ are cosine of the reconstructed 
zenith angle of muon and neutrino, respectively. " 9 J 0- 

In our understanding, SK Monte Carlo Simulation is 
named usually as the Detector Simulation. It is, however, 
noticed that the effect of the azimuthal angles of the emit- 
ted leptons in QEL could not be taken into account in their 
Simulation. As will be shown in later (see Section 3), this 
effect greatly influences over the final zenith angle distri- 
bution of the emitted leptons. Also, the back scattering 
due to QEL can not be neglected for the rigorous deter- 
mination of the direction of the incident neutrino, but this 
effect could not be treated in the SK Detector Simulation, 
which is beyond the application limitaition 0. 

On the other hand, we could take into account these 
effects correctly in our Monte Carlo Simulation which is 
named as Time Seguential Simulation. 

In the present paper, we carry out the full Monte Carlo 
Time Sequential Simulation as exactly as possible, with- 
out the SK assumption on the direction to clarify the prob- 
lematic issue raised by SK. We carry out simulation which 
starts from the opposite side of the Earth to the SK de- 
tector. A neutrino sampled from the atmospheric neutrino 
energy spectrum at the opposite side of the Earth tra- 
verses through the medium with different densities in the 
interior of the Earth and penetrates finally into the SK 



5 It should be noticed that the SK assumption on the di- 
rection may hold on the following possible two cases: [1] The 
scattering angle of the emitted lepton is so small that the ef- 
fect of the scattering angle could be neglected really. However, 
in the present case, it could not be true from Fig. [T]and Fig. 
[2] and Table 1. [2] One may assert that the assumption could 
not hold on individual case, but it could hold statistically af- 
ter accumulation of large amount of the data. However, such 
assertion should be verified. We verify such assumption could 
not hold. See, Fig. 1111 and Fig. [12] in the present paper. 

6 SK have never clarified not only the details, but also the 
principle and its validity on their Monte Carlo Simulation. We 
hope disclosure of their Detector Simulation for open and fair 
scientific discussion. 



detector where the neutrino interactions occur. The emit- 
ted energy of the individual lepton thus produced and its 
direction are simulated exactly based on the probability 
function of the cross sections concerned. 

We finally show the zenith angle distribution of the 
emitted leptons as well as that of the incident neutrinos 
are quite different from corresponding ones of the SK. This 
indicates that the SK assumption on the direction coud 
not be a reliable estimator as for the determination of the 
direction of the incident neutrino (See, section 5). 



2 Cross Sections of Quasi Elastic Scattering 
in the Neutrino Reaction and the Scattering 
Angle of Charged Leptons. 

We examine the following reactions due to the charged 
current interaction (c.c.) from QEL. 
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The differential cross section for QEL is given as fol- 
lows (6l. 
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The signs + and — refer to ^(e) and ^W e ) for charged 
current (c.c.) interactions, respectively. The Q 2 denotes 
the four momentum transfer between the incident neu- 
trino and the charged lepton. Details of other symbols are 
given in [4j. 

The relation among Q 2 , E v r p \, the energy of the in- 
cident neutrino, Ej>, the energy of the emitted charged 
lepton (muon or electron or their anti-particles) and 9 S , 
the scattering angle of the emitted lepton, is given as 



Q 2 = 2£ v(p) £k(l-cos0 s 



(3) 
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Fig. 1. Relation between the energy of the muon and its 
scattering angle for different incident muon neutrino energies, 
0.5 GeV, 1 GeV, 2 GeV, 5 GeV, 10 GeV and 100 GeV. 



Also, the energy of the emitted lepton is given by 




Fig. 2. Distribution functions for the scattering angle of the 
muon for muon-neutrino with incident energies, 0.5 GeV, 1.0 
GeV and 2 GeV. Each curve is obtained by the Monte Carlo 
method (one million sampling per each curve). 



Now, let us examine the magnitude of the scattering 
angle of the emitted lepton in a quantitative way, as this 
plays a decisive role in determining the accuracy of the 
direction of the incident neutrino, which is directly related 
to the reliability of the zenith angle distribution of both 
Fully Contained Events and Partially Contained Events in 
SK. 

By using Eqs. ([2|) to (JI]), we obtain the distribution 
function for the scattering angle of the emitted leptons 
and the related quantities by a Monte Carlo method. The 
procedure for determining the scattering angle for a given 
energy of the incident neutrino is described in the Ap- 
pendix A. Fig. [T] shows this relation for muon, from which 
we can easily understand that the scattering angle 9 S of 
the emitted lepton ( muon here ) cannot be neglected. 
For a quantitative examination of the scattering angle, 
we construct the distribution function for 6 S of the emit- 
ted lepton from Eqs. ([2]) to ([4| by using a Monte Carlo 
method. 

Fig.[2]gives the distribution function for 9 S of the muon 
produced in the muon neutrino interaction. It can be seen 
that the muons produced from lower energy neutrinos are 
scattered over wider angles and that a considerable part 
of them are scattered even in backward directions. Simi- 
lar results are obtained for anti-muon neutrinos, electron 
neutrinos and anti-electron neutrinos. 

Also, in a similar manner, we obtain not only the dis- 
tribution function for the scattering angle of the charged 
leptons, but also their average values < 9 S > and their 
standard deviations o~ s . Table 1 shows them for muon neu- 
trinos, anti-muon neutrinos, electron neutrinos and anti- 
electron neutrinos. In the SK analysis, it is assumed that 
the scattering angle of the charged particle is zero [8j£9] . 



3 Influence of Azimuthal Angle of Quasi 
Elastic Scattering over the Zenith Angle of 
both the Fully Contained Events and 
Partially Contained Events 



In the present section, we examine the effect of the az- 
imuthal angles of the emitted leptons over their own zenith 
angles for given zenith angles of the incident neutrinos Q. 

For three typical cases (vertical, horizontal and diag- 
onal), Fig. 3 gives a schematic representation of the re- 
lationship between, #„(p), the zenith angle of the incident 
neutrino, and (9 S , </>) a pair of scattering angle of the emit- 
ted lepton and its azimutal angle. 

From Fig. 3(a), it can been seen that the zenith angle 
of the emitted lepton is not influenced by its <j) in 
the vertical incidence of the neutrinos {0 v (v) — 0°), as it 
must be. From Fig. 3(b), however, it is obvious that the 
influence of <fi of the emitted leptons on their own zenith 
angle is the strongest in the case of horizontal incidence 
of the neutrino {9 v (v) = 90°). Namely, one half of the 
emitted leptons are recognized as upward going, while the 
other half is classified as downward going ones. The di- 
agonal case ( v (v) — 43°) is intermediate between the 
vertical and the horizontal. In the following, we examine 
the cases for vertical, horizontal and diagonal incidence of 
the neutrino with different energies, say, E v rp\ — 0.5 GeV, 
E v{p) = 1 GeV and E v{p) = 5 GeV. 

The detailed procedure for the Monte Carlo simulation 
is described in the Appendix A. 



7 Throughout this paper, we measure the zenith angles of 
the emitted leptons from the upward vertical direction of the 
incident neutrino. Consequently, notice that the sign of our 
direction is oposite to that of the SK ( our cos Q v (p) = - cos Q v (v) 
in SK) 
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Table 1. The average values < 9 S > for scattering angle of the emitted charged leptons and their standard deviations o s for 
various primary neutrino energies E v ^. 
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Fig. 3. Schematic view of the zenith angles of the charged 
muons for diffrent zenith angles of the incident neutrinos, fo- 
cusing on their azimuthal angles. 



3.1 Dependence of the spreads of the zenith angle for 
the emitted leptons on the energies of the emitted 
leptons for different incident directions with different 
energies 



We give the scatter plots between the fractional energies 
of the emitted muons and their zenith angle for a definite 
zenith angles of the incident neutrino with different ener- 
gies in Figs. [3] to In Fig. 0] we give the scatter plots 
for vertically incident neutrino with different energies 0.5 
GeV, 1 GeV and 5 GeV . In this case, the relations between 
the emitted energies of the muon and and their zenith an- 
gles are unique, which comes from the definition of the 
zenith angle of the emitted lepton. However, the densities 
(frequencies of event number) along each curve is differ- 
ent in position to position and depend on the energies of 
the incident neutrinos. Generally speaking, densities along 
curves become higher toward cos = 1. In this case, 
cos O^tjA is never influenced by the azimuthal angel in the 
scattering by the definition 0. 

Fig. 5 tells us that the horizontally incident neutrinos 
give the most widely spread of the zenith angle distribu- 
tion of the emitted lepton influenced by the azimuthal an- 
gle. The more lower incident neutrino energies, the more 
wider spreads of the emitted leptons. The diagonally in- 
cident neutrinos give the intermediate distribution of the 
emitted leptons between those of vertically incident neu- 
trinos and horizontally incident neutrinos. 



The zenith angles of the particles concerned are measured 
from the vertical direction. 



6 E.Konishi et. al.,: The Reliability on the Direction of Neutrino in SK 




E,/E v E,/E v E^/E v 

Fig. 4. The scatter plots between the fractional energies of the produced muons and their zenith angles for vertically incident 
muon neutrinos with 0.5 GeV, 1 GeV and 5 GeV, respectively. The sampling number is 1000 for each case. 
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Fig. 5. The scatter plots between the fractional energies of the produced muons and their zenith angles for horizontally incident 
muon neutrinos with 0.5 GeV, 1 GeV and 5 GeV, respectively. The sampling number is 1000 for each case. 
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Fig. 6. The scatter plots between the fractional energies of the produced muons and their zenith angles for diagonally incident 
muon neutrinos with 0.5 GeV, 1 GeV and 5 GeV, respectively. The sampling number is 1000 for each case. 
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Fig. 7. Zenith angle distribution of the muon for the vertically incident muon neutrino with 0.5 GeV, 1 GeV and 5 GeV, 
respectively. The sampling number is 10000 for each case. SK stand for the corresponding ones under the SK assumption. 
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Fig. 8. Zenith angle distribution of the muon for the horizontally incident muon neutrino with 0.5 GeV, 1 GeV and 5 GeV, 
respectively. The sampling number is 10000 for each case. SK stand for the corresponding ones under the SK assumption. 



(a) 



mu 

CO! 


Dn r 

;9 V = 


eut 
0.7, 


'ino 
31 (6 


v=4 


3°) 










a\ 


'era 


OrW 
ge= 


3e\ 
0.1 


39 






s 


K 






s 


d.= 


0.5£ 


56 





























































































-0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 
COS9„ 



CD 
CO 

o 
o 

"D 



(b) 



mu 

CO! 


Dn r 

;9 V = 


eut 
0.7, 


'i no 
31(6 


v=4 


3°) 










a\ 


'era 


H3i 
ge= 


3V~ 

0.4: 


32 






s 


K 






s 


d.= 


0.46 


S3 









































































-0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 
COS9„ 



CD 
CO 

o 
o 



muon neut 
co?9 v 



average 



(c) 

ino 

0.7Sl(9| v =4f3 c 
5GeV 



0.715 



0.103 



-0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 
COS9„ 



Fig. 9. Zenith angle distribution of the muon for the diagonally incident muon neutrino with 0.5 GeV, 1 GeV and 5 GeV, 
respectively. The sampling number is 10000 for each case. SK stand for the corresponding ones under the SK assumption. 
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3.2 Zenith angle distribution of the emitted lepton for 
the different incidence of the neutrinos with different 
energies 

In Figs. [7] to O we give the zenith angle distributions of 
the emitted muons for the given direction of the incident 
neutrinos with different energies of the neutrino. These 
figures are obtained through summation on the energies 
of the emitted muons for their definite zenith angles in 
Figs. U to E 

In Figs. [7ja) to [2c), we give the zenith angle distri- 
bution of the emitted muon for the case of vertically inci- 
dent neutrinos with different energies, say, E v — 0.5 GeV, 
E v = 1 GeV and E v = 5 GeV. 

Comparing the case for 0.5 GeV with that for 5 GeV, 
we understand the big contrast between them as for the 
zenith angle distribution. The scattering angle of the emit- 
ted muon for 5 GeV neutrino is relatively small (See, Table 
1) that the emitted muons keep roughly the same direction 
as their original neutrino. In this case, the effect of their 
azimuthal angle on the zenith angle is also small. However, 
in the case for 0.5 GeV which is the dominant energy for 
Fully Contained Events in the Superkamiokande, there is 
even a possibility for the emitted muon to be emitted in 
the backward direction due to the large angle scattering, 
the effect of which is enhanced by their azimuthal angle. 

The most frequent occurrence in the backward scatter- 
ing of the emitted muon appear in the horizontally inci- 
dent neutrino as shown in Figs. 8(a) to 8(c). In this case, 
the zenith angle distribution of the emitted muon should 
be symmetrical to the horizontal direction. Comparing the 
case for 5 GeV with those both for ~0.5 GeV and for ~1 
GeV, even 1 GeV incident neutrinos lose almost the orig- 
inal sense of the incidence if we measure it by the zenith 
angle of the emitted muon. Figs. 9(a) to 9(c) for the di- 
agonally incident neutrino tell us that the situation for 
diagonal cases lies between the case for the vertically in- 
cident neutrino and that for horizontally incident one. 

4 Zenith Angle Distribution of Fully 
Contained Events and Partially Contained 
Events for a Given Zenith Angle of the 
Incident Neutrino, Taking Their Energy 
Spectrum into Account 

In the previous sections, we discuss the relation between 
the zenith angle distribution of the incident neutrino with 
a single energy and that of the emited muons produced by 
the neutrino for the different incident direction. In order to 
apply our motivation around the uncertainty of the SK as- 
sumption on the direction for Fully Contained Events and 
Partially Contained Events, we must consider the effect of 
the energy spectrum of the incident neutrino. The Monte 
Carlo simulation procedure for this purpose are given in 
the Appendix B. 

In Fig. [TD1 we give the zenith angle distributions of the 
sum of jJL + (fi) and (JT for a given zenith angle of and 



z^j, taking into account primary neutrino energy spectrum 
at Kamioka site. 

In Table 2, the average values for cos6> M+ ^ and their 
standard deviation for different incidences of the incident 
neutrinos with different energies are presented 0- 

In the SK case, their average values are given by coa6 v ^ 
themselves by definition and, consequently, the standard 
deviations are zero under the assumption, because the SK 
assumption is of the delta function for the incidence direc- 
tion. They are shown in the bottom line of Table 2. In the 
second line from the bottom in this table, we give the av- 
erage values and their standard deviations for cos#„+ p ob- 
tained under the inclusion of the energy spectrum for pri- 
mary neutrinos. Thus, we found these values correspond to 
those for incident neutrino with the effective single energy 
between 0.5 GeV and 1 GeV. If we compare the average 
energies and the standard deviations for the inclusion of 
incident neutrino energy spectrum with those under the 
SK assumption, it is easily understood that SK assump- 
tion does not represent real zenith angle distribution of 
the emitted muon. 

5 Relation between the Zenth angle 
Distribution of the Incident Neutrinos and 
that of the emitted leptons 

Now, we extend the results for the definite zenith angle 
obtained in the previous section to the case in which we 
consider the zenih angle distribution of the incident neu- 
trinos totally. 

Here, we examine the real correlation between cos 6 V 
and cos# M , by peforming the exact Monte Carlo simula- 
tion. 

The detail for the simulation procedure is given the 
Appendix C. 

In Fig. [TTJ we classsify the correlation between cos B v 
and cos 6^ according to the different energy range of the 
incident muon neutrinos. It should be noticed that the 
SK assumption on cos B v = cos 0^ is roughly hold only for 
E v > 5 GeV, but the widths in cos 9^ for the definite cos 6 V 
near cos 6 V =0 (for horizontally incident neutrino ) are 
much larger than those near cosf^ =1 (for the vertically 
incident neutrino). Of course, this is due to the effect of 
the azimutal angle in QEL which could not be derived by 
the SK simulation (DETECTOR SIMULATION). Such 
tendencies become more remarkable in E v < 5 GeV and 
in these energies the SK assumption on the direction does 
not hold any more. 

In Fig. Q21 we classify the correlation between cos d v and 
cos Ofj, according to the different energy range of . The 
similar argument on Fig. [TT] can be done on the case of 
FiglU 

9 Notice that the difference in the corresponding quantites 
between the case for single energy and the case for the energy 
spectrum. The formers are given in the u~ , while th latter is 
given in u~ and u + . However, such the difference does not 
change the essential recognition. 
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Table 2. Average values and their standerd deviations in cosO^+p for the zenith angle distributions of the muons with different 
primary energies of the insident neutrinos. 
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Fig. 10. Zenith angle distribution of and /i + for v and v for the incident neutrinos with the vertical, horizontal and diagonal 
directions, respectively. The overall neutrino spectrum at Kamioka site is taken into account. The sampling number is 10000 
for each case. SK stand for the corresponding ones under the SK assumption. 
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Thus, it could be surely concluded from Fig. QT] and 
Fig. [T2] that the SK assumption on the direction never 
holds as a good estimator for the determination of the 
directions of the incident neutrinos. 

In order to obtain the zenith angle distribution of the 
emitted leptons for that of the incident neutrinos, we di- 
vide the cosine of the zenith angle of the incident neutrino 
into twenty regular intervals from cos 9 V = to cos 6 V = 1 . 
For the given interval of cos 6 U , we carry out the exact 
Monte Carlo simulation, the detail of which is give in the 
Appendix D and obtain the cosine of the zenith angle of 
the emitted leptons, taking account of the geometry for 
surronding the SK detector. 

Thus, for each interval of cos^, we obtain the corre- 
sponding zenith angle distribution of the emitted leptons. 
Then, we sum up these corresponding ones over all zcnth 
angles of the incident neutrinos and we finally obtain the 
relation between the zenith angle distribution for the in- 
cident neutrinos and that for the emitted leptons. 

In a similar manner, we could obtain between cos 9p 
and cos 9 fi for anti-neutrinos. The situation for anti-neutrinos 
is essentially same as that for neutrinos. 

Here, we examine the zenith angle distribution of the 
muons from both upward neutrinos and downward ones 
in the case that neutrino oscillation does not exist. 

By performing the procedures described in Appendix 
C, a pair of sampling ( cos^+p, E w+j y ) gives a pair of ( 
cos 9fi+fi , Ep+p. ) ■ In Fig. [13l we give the zenith angle dis- 
tribution of the upward neutrinos ( the sum of and Pp ) 
which is constructed from the energy spectra for different 
cos#j, +p . (see, Honda et. al. fTU] and Appendix B) 

Upward neutrinos may produce even downward lep- 
tons due to both the backscattering effect and the effect 
of azimuthal angle on larger forward scattering for the in- 
teraction concerned (see, Figure 3 and Figures 4 to 6 in 



the text). As the result of it, the zenith angle distribution 
of the emitted muons for the upward neutrino may leak in 
the downward direction. From Figure 13, it is very clear 
that the shape of the zenith angle distribution for the in- 
cident neutrinos is quite different from that of the emitted 
muons produced by these neutrinos. If the SK assumption 
on the direction statistically holds, the zenith angle dis- 
tribution for the emitted muons should coincide totally 
with that of the incident neutrinos. In other words, one 
may say that the zenith angle distribution for the emitted 
muons should be understood as that of the incident neu- 
trino under the SK assumption on the direction. However, 
the muon spectrum is distinctively different from the real 
(computational) incident neutrino spectrum as shown in 
the figure. Thus we conclude that SK assumption on the 
direction leads to the wrong conclusion on the neutrino 
oscillation. The further examination on the experimental 
data obtained by SK will be carried out in the subsequent 
papers. 

It is, further, noticed that upward neutrino energy 
spectrum in the figure biggest near cos = and the 
smallest near cos(9„ +p = 1, which reflects from the en- 
hancement of the primary incident neutrino energy spec- 
trum from the inclined direction and is independent on the 
neutrino oscillation, while in SK opinion, such tendency 
may be favor of the existence of neutrino oscillation. 



6 Discussions and Conclusion 

In order to extract the definite conclusion on the neutrino 
oscillation from the experiment by cosmic ray neutrinos 
whose intensity as well as interaction with the substance 
are both very weak, first of all, one should analyze the 
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most clear cut and ambiguity free events. Among neutrino 
events analyzed by SK, the most clear cut events are Single 
Ring Events, such as electron-like events and muon-like 
events in Fully Contained Events which are generated by 
QEL. These events are provided with simplicity due to 
single ring and all possible measureable physical quantities 
are confined in the detector. 

Furthermore, QEL is the most dominant source for 
neutrino events which are generated in the SK detector. 
This is the reason why we examine the QEL events excul- 
sively in present and subsequent papers. 

If the neutrino oscillation really exists, the most clear 
cut evidence surely appears in the analysis of single ring 
events due to QEL in Fully Contained Events and one does 
not need the analysis of any other type of events, such as 
single ring events in Partially Contained Events, multi- 
ring events in either Fully Contained Events or Partially 
Contained Evetns, all of which include inevitably ambigu- 
ities for the interpretation and show merely sub-evidences 
compared with that from the single ring events due to 
QEL in Fully Contained Events. 

SK analyze the zenith angle distribution of the inci- 
dent neutrinos under the asssumption that the direction 
of the incident neutrino is the same as that of the emitted 
lepton. We conclude that this assumption is supplemented 
by their Monte Carlo Simulation named as Detector Sim- 
ulation^^. 

In the present paper, we adopt Time Sequential Simu- 
lation which starts from the incident neutrino energy spec- 
trum on the opposite side of the Earth to the SK detec- 
tor and simulate all posible physical processes which are 
connected with the zenith angle distribution of the in- 
cident neutrinos according to their probability functions 
concerned for the examination on the validity of the SK 
assumption on the direction. 

Concretely speaking, we take the following treatment, 
(i) the stochastic treatments of the scattering angle of the 
emitted lepton in QEL, including the scattering on the 
backward as well as the azimuthal angle, which could not 
be treated in Detector Simulation, (ii) the stochastic treat- 
ment on the zenith angle distribution of the emitted lep- 
ton, considering the incident neutrino energy spectrum, 
(hi) the stochastic treatment on the detection of the QEL 
events inside the SK detector. 

Furthermore, the discrimination between Fully Con- 
tained Events and Partially Contained Events is only pos- 
sible in the Time Sequential Simulation, because the events 
concerned may be classified into different categories by 
chance, Fully Contained Events and Partially Contained 
Events due to different occurring points and different di- 
rections. 



The SK Detector Simulation for obtaing the zenith angle 
distribtuion of the incident neutrino as for the neutrino oscilla- 
tiom has never been disclosed in their papers, even in the Ph.D 
thesis. Consequently, this is only our onjecture as for utiliza- 
tion of the SK Detector Simulation . A clear thing is only that 
SK impose the proposittion that the direction of the incident 
neutrino is the same as that of the emitted lepton upon the 
neutrino oscillation analysis. 



The conclusions thus obtained are as follows: 

(1) The zenith angle distributions of the emitted lepton 
in QEL for the incident neutrino with both the def- 
inite zenith angle and the definite energy are widely 
spread, particularly, into even the backward region due 
to partly pure backscattering and partly the combina- 
tion of the azimuthal angle with the slant direction 
of the incident neutrinos. However, for every incident 
neutrino with a definte zenith angle, SK give the same 
definite zenith angle to the emitted lepton. Already in 
this stage, the SK assumption on the direction does 
not hold. 

(2) Taking account of the incident neutrino energy spec- 
trum and simulating all physical processes concerned, 
we obtain the zenith angle distribution of the emitted 
leptons for the incident neutrino with a definite zenith 
angle. It is proved that the SK assumption on the di- 
rection does not hold again. 

(3) The correlation diagrams between cosd u and cos6> M 
show that SK assumption does noty hold well even 
for higher energies of the incident neutrinos, and it 
is shown that the correlation between them become 
weaker in more inclined incident neutrinos due to the 
effect of the azimuthal angle in QEL. 

(4) Taking into account the detection efficiency for the 
events concerned in the simulation for upward neu- 
trinos and anti-neutrinos, we obtain the zenith angle 
distribution of the leptons ( muons plus anti-muons 
). According to the SK Assumption on the direction, 
the zenith angle distribution is the same as that of the 
incident neutrinos. However, the original zenith angle 
distribution of incident neutrino is found to be quite 
different from that derived from that of leptons. This is 
the final conclusion that SK have not measured the di- 
rection of the incident neutrinos reliably, which is quite 
independent on either the existence or non-existence of 
the neutrino oscillation. 

(5) The SK assume that the Partially Contained Events 
exclusively belong to the muon-like event. However, 
such the assumption lacks in theoretical background. 
Electron events can also contribute to the Partially 
Contained Events under some geometrical condition, 
for example, partly coming from the transformation 
by Eq.(A.5). The quantitative examination on the Par- 
tially Contained Events among the electron-like event 
will be published elsewhere. 

In subsequent papers, we will give the relation between 
the zenith angle distributions of the incident neutrinos 
and the corresponding muons in the cases with and with- 
out neutrino oscillation, including downward neutrino and 
will examine whether it is possible to or not to detect the 
neutrino oscillation by using atmospheric neutrino. 



12 



E.Konishi et. al.,: The Reliability on the Direction of Neutrino in SK 



In the following Appendices we give the concrete Monte 
Carlo Simulations, namely, the details of our Time Se- 
quential Simulation. 



A Appendix: Monte Carlo Procedure for the 
Decision of Emitted Energies of the Leptons 
and Their Direction Cosines 

Here, we give the Monte Carlo Simulation procedure for 
obtaining the energy and its direction cosines, (l r , m r , n r ), 
of the emitted lepton in QEL for a given energy and its 
direction cosines, (l,m,n), of the incident neutrino. 

The relation among Q 2 , E u +p, the energy of the inci- 
dent neutrino, Et, the energy of the emitted lepton (muon 
or electron or their anti-particles) and 9 S , the scattering 
angle of the emitted lepton, is given as 

Q 2 = 2E v{v) E m (\ - cosfl s ). (AT) 

Also, the energy of the emitted lepton is given by 



9L 

2M' 



(A-2) 



Procedure 1 

We decide Q 2 from the probability function for the differ- 
ential cross section with a given E v i„\ (Eq. (|2|) in the text) 
by using the uniform random number, £, between (0,1) in 
the following 



Q 



Pt(i)(Ev(„),Q 2 )dQ 2 



(A-3) 



where 

Pl(l){ E „(y),Q ) = 

&Vl{Z){E u {p)i Q 2 ) 
dQ 2 




Q,iax &<j l{l) {E v{i>) ,Q 2 ) 



dQ 2 



dQ'. 
(A-4) 




0.4 0.8 

Q 2 /Q 



max 



Fig. 14. The reappearance of the probability function for QEL 
cross section. Histograms are sampling results, while the curves 
concerned are theoretical ones for given incident energies. 



which is obtained from 



= 2<. 



(A-5) 



Here, £ is a uniform random number (0, 1). 
As explained schematically in the text (see Fig. [3] in the 
text), we must take account of the effect due to the az- 
imuthal angle 4> in the QEL to obtain the zenith angle 
distribution of both Fully Contained Events and Partially 
Contained Events correctly. 

Procedure 5 

The relation between direction cosines of the incident neu- 
trinos, (4(j/)! m i/(j/)i n i/(i/))i an< i those of the corresponding 
emitted lepton, (£ r , m x , rt r ), for a certain 8 S and <f> is given 
as 



From Eq. (A-l), we obtain Q 2 in histograms together with 
the corresponding theoretical curve in Fig. Q3] The agree- 
ment between the sampling data and the theoretical curve 
is excellent, which shows the validity of the utlized proce- 
dure in Eq. (A-3) is right. 

Procedure 2 

We obtain E t ^ from Eq. (A-2) for the given E v ^ and 
Q 2 thus decided in the Procedure 1. 

Procedure 3 

We obtain cos 6> s , cosine of the the scattering angle of the 
emitted lepton, for Etg\ thus decided in the Procedure 2 
from Eq. (A-l) . 

Procedure 4 

We decide cj>, the azimuthal angle of the scattering lepton, 



in 









- 







\Jl 2 + m 2 



VP 



VP 



(sin^scost/; 
sin6* s sin(/) 
cos6> s , 

(A-6) 



where n 



cos8 



and n r = cos6*^. Here, 0£ is the 



zenith angle of the emitted lepton. 

The Monte Carlo procedure for the determination of Bg 
of the emitted lepton for the parent (anti-) neutrino with 
given Q v (jj) and £v( P ) involves the following steps: 



We obtain (£ r ,m ri n r ) by using Eq. (|A-6[) . The n r is 
the cosine of the zenith angle of the emitted lepton which 
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Fig. 15. The relation between the direction cosine of the 
incident neutrino and that of the emitted charged lepton. 



should be contrasted to n v , that of the incident neutrino. 
Repeating the procedures 1 to 5 just mentioned above, we 
obtain the zenith angle distribution of the emitted leptons 
for a given zenth angle of the incident neutrino with a def- 
inite energy. 



In the SK analysis, instead of Eq. (|A 6|) . they assume 
n r = n u i V ) uniquely for > 400 MeV. 



B Appendix: Monte Carlo Procedure to 
Obtain the Zenith Angle of the Emitted 
Lepton for a Given Zentith Angle of the 
Incident Neutrino 

The present simulation procedure for a given zenith an- 
gle of the incident neutrino starts from the atmospheric 
neutrino spectrum at the opposite site of the Earth to 
the SK detector. We define, Ni a t{E u i v \,t, cos#„( P )), the in- 
teraction neutrino spectrum at the depth t from the SK 
detector in the following way 

^Vint(-E ; I /(p),*,COS6' 1 ,(p)) = N sp (-E„(p), COS 0„(p))X 

dt \ ( dt \ 



1 - 



AiX-Epfp), h, pi) 



x • • • x 1 - 



Xn{E 

u{u) 1 ^"ni Pn) J 

(B-l) 



Here, N sv {E u ^ i} ) 1 cos9 u ^)) is the atmospheric (anti-) 
neutrino spectrum for the zenith angle at the opposite 
surface of the Earth. 

Here Xi(E v ^,ti, pi) denotes the mean free path due 
to the neutrino(anti neutrino) with the energy E V ( P \ from 
QEL at the distance, U, from the opposite surface of the 
Earth inside whose density is pi. 

The procedures of the Monte Carlo Simulation for the 
incident neutrino(anti neutrino) with a given energy, E v ip\ , 
whose incident direction is expressde by (I, m, n) are as fol- 
lows. 

Procedure A 

For the given zenith angle of the incident neutrino, 9 U ^ 1 
we formulate, N vm {E v ^, i, cos9 v ^)dE v ^, the produc- 
tion function for the neutrino flux to produce leptons at 
the Kamioka site in the following 

-^pro(-^V(p) i ^5 COS 9 v (p} 

= (T i(l){Ev(v))Nint{E v ( U - )1 t,cos9 v ^)dE u ^i (B-2) 



where 

a e(e){E u (p)) 



^t{l){Eu(v)iQ 2 ) 
dQ2 



dQ 2 



(B-3) 



Each differential cross section above is given in Eq. ((2|) in 
the text. 

Utilizing, £, the uniform random number between (0,1), 
we determine -E„( P ), the energy of the incident neutrino in 
the following sampling procedure 



where 



Pd{E u{p) , t, cos 0„ ( p) (v))dE v{p) , (B-4) 



Pd{Ev(p} , t, cos 9 v (p^ )d^ ( P) 

N P ro{E v ^) , t, cos 6>^ (P ) )dE v ^) 

.(ff),max 

Npro{E^^),t, COS ^y(p))di? l/ (p) 



(B-5) 
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In our Monte Carlo procedure, 
the reproduction of, Pd(E v tp\,t, cos d^tn^dE^f^, the nor- 
malized differential neutrino interaction probability func- 
tion, is confirmed in the same way as in Eq. (A-4). 



Procedure B 

For the (anti-)neutrino concerned with the energy of -E„( P ) , 
we sample Q 2 utlizing £3, the uniform random number be- 
tween (0,1). The Procedure B is exactly the same as in the 
Procedure 1 in the Appendix A. 

Procedure C 

We decide, 9 S , the scattering angle of the emitted lepton 
for given E v i v \ and Q 2 . The procedure C is exactly the 
same as in the combination of Procedures 2 and 3 in the 
Appendix A. 

Procedure D 

We randomly sample the azimuthal angle of the charged 
lepton concerned. The Procedure D is exactly the same as 
in the Procedure 4 in the Appendix A. 

Procedure E 

We decide the direction cosine of the charged lepton con- 
cerned. The Procedure E is exactly the same as in the 
Procedure 5 in the Appendix A. 

We repeat Procedures A to E until we reach the de- 
sired trial number. 



C Appendix: Correlation between the Zenith 
Angles of the Incident Neutrinos and Those 
of the Emitted Leptons 

Procedure A 

By usi ng, N plo (E l/ ( j? ) , t, cos B v (^ )dE l/(P ) , which is defined 
inEq. (ED, 

we define the spectrum for cos in the following. 
7(cos6' l/ (p))d(cos^(p ) ) = 

max 

d(cos6» y ( P )) / N VTO (E v (p),t, cos6» iy( p))d£; iy (p).(C-l) 
JE uiP)>tDia 

By using Eq. (|C-2[) and £ , a sampled uniform random num- 
ber between (0,1), then we could determine cos#„(p) from 
the following equation 

P n (cos 6» l/(p) )d(cos d v {p)), (C-2) 



/■COS t 



where 

P n (cos 0^(5)) = /(cos v (v)) j J /(cos6» iy(P ))d(cos^(p)). 

(C3) 

Procedure B 

1-1 ii) n j) in the Procedure A, we sam- 

using £, the uniform randum 



Procedure B 

For the sampled 6\{cosd v (p\) in the Procedure A, wc 
pie E v fp\ from Eq. (IC-4l) by using £, the uniform ra 
number between (0,1) 



Ppro{E v (v) , COS Q v {p} )d-E„(p) , (C-4) 



where 



Ppro{E v (v),t,COs9v{v))&E v ( 9 ) — 
•^pro(-^V(p) ) t : COS u ^v) 

)dE v 
Procedure C 

For the sampled E v ^ in the Procedure B, we sample 
Eptft from Eqs. (|A-2I) and (|A-3[) . For the sampled E v i 9 \ 
and E^(n\, we determine cos0 s , the scattering angle of the 
muon uniquely from Eq. (|A-1[) . 

Procedure D 

Wc determine, <fi, the azimuthal angle of the scattering lep- 
ton from Eq. (|A-5|) by using £, an uniform randum number 
between (0,1). 

Procedure E 

We obtain cos from Eq. (|A-6|) . As the result, we ob- 
tain a pair of (cos cos# M (p)) through Procedures A 
to E. Repeating the Procedures A to E, we finally the cor- 
relation between the zenith angle of the incident neutrino 
and that of the emitted muon. 
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D Appendix: Detection of the Neutrino 
Events in the SK Detector and Their 
Interaction Points 

The plane ABCD is always directed vertically to the di- 
rection of the incident neutrino with a given zenith an- 
gle, which is shown in Fig. [16] The rectangular ABCDE- 
FGH encloses the SK detector whose radius and height 
is denoted by R and H, respectively. The width and the 
height of the plane ABCD for a given zenith angle, 0„(p) , 
is given as, R and Rcos9 u ^ + Hsva9 v , respectively, which 
are shown in Fig. 16-c. 

Now, let us estimate the ratio of the number of the neu- 
trino events inside the SK detector to that in the rectan- 
gular ABCDEFGH. As the number of the neutrino events 
inside some material is proportional to the number of the 
nucleons in the material concerned. The number of the 
nucleons inside the SK detector (p = 1) is given as 



./VextrtcOS^p) ) = pNg, 



1 - 



R 2 H 



R(H 2 +R 2 ) sin 9 V ^ cos d v (p) 



, (D-2) 



where p is the density of the rock which surrounds the SK 
detector. 

Then, the total number of the target in the rectangu- 
lar ABCDEFGH is given as 



N toi (cos6 v (y)) = N sk + iV extr (cos0„(p) 
Here, we take 2.65, as p (standard rock). 



(D-3) 



Then, i?thoor, the ratio of the number of the neutrino 
events in the SK detector to that in the rectangular ABCDE- 
FGH is given as 



N sk = jN avoga R 2 H, (D'l) 

where N avoga denotes the Avogadro number, and the num- 
ber of the nucleons in the exterior of the SK detector inside 
ABCDEFGH is given as 




^thcor 

(cos#„( P )) = iV s k/-/Vtot (cos ^(p)). (D-4) 

We obtain -Rtheor for different values of cos^p) given in 
the Table 3. 

Here, we simulate neutrino events occured in the rectangu- 
lar ABCDEFH, by using the atmospheric neutrino beam 
which falls down on the plane ABCD. Thus, iV sma pi e , the 
sampling number of the (anti-)neutrino events inside the 
rectangular ABCDEFG for a given cos6 u ^ is given as 

^sample (COS v{p )) = N to t (COS v{9) ) X 

/ a M) (E v (j?))N int {E v{p ) , t, cosfl^p) )&E v {p) 

(D-5) 

where ^ (E^^) is the total cross section for (anti-)neutrino 
due to QEL, and N- m t{E v ^,t,co86 v [p))dE v {p) is the dif- 
ferential nutrino energy spectrum for the definite zenith 
angle, 9 v fp), in the plane ABCD. The injection points of 
the neutrinos in the plane ABCD are distributed over the 
plane randomly and uniformely and the injection points 
are determined from a pair of the uniform random num- 
bers between (0,1). They penetrate into the rectangular 
ABCDEFGH from the injection point in the plane ABCD 
and some of them may penetrate into the SK detector or 
may not, which depend on their injection point. 

In the neutrino events which penetrate into the SK 
detectorr, their geometrical total track length, T trac k, are 
devided into three parts 

Ttrack = T h + T sk + T a , (D-6) 



Fig. 16. Sampling procedure for neutrino evens injected into 
the detector 



where Tb denotes the track length from the plane ABCD 
to the entrance point of the SK detector, T s k denotes the 
track length inside the SK detector, and T a denotes the 
track length from the escaping point of the SK detector 
to the exit point of the rectangular ABCDEF, and thus 
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Table 3. Occurrence probabilities of the neutrino events in- 
side the SK detector for different cos 9 V 's. Comparison between 
■Rthcor and -Rmontc- The sampling numbers for the Monte Carlo 
Simulation are, 1000, 10000, 100000, respectively. 



cos f5„ 



R 



th eor 



Sampling Number 



1000 10000 100000 



0.000 





58002 





576 





5750 


0.57979 


0.100 





41717 





425 





4185 


0.41742 


0.200 





32792 





353 





3252 


0.32657 


0.300 





27324 





282 





2731 


0.27163 


0.400 





23778 





223 





2329 


0.23582 


0.500 





21491 





206 





2063 


0.21203 


0.600 





20117 





197 





1946 


0.19882 


0.700 





19587 





193 





1925 


0.19428 


0.800 





20117 





198 





2002 


0.20001 


0.900 





22843 





230 





2248 


0.22803 


1.00 





58002 





557 





5744 


0.57936 



Ttrack denotes the geometrical length of the neutrino con- 
cerned in the rectangular ABCDEFGH. 

By the definition, the neutrinos concerned with T trac k 
interact surely somewhere along the Ttrack- Here, we are 
interested only in the interaction point which ocuurs along 
T sk . We could determine the interaction point in the T sk 
in the following. 

We define the following quantities for the purpose. 



weight 
Pav 



= T sk + p(T b + T a ), 



T 

1 w 



Ait 



/T tr 



acki 



= Pav/P, 



2-ht- 



(D-7) 
(D-8) 
(D-9) 
(D-10) 



The flow chart for the choice of the neutrino events in 
the SK detector and the determination of the interaction 
points inside the SK detector is given in Fig. [171 Thus, we 
obtain neutrino events whose occurrence point is decided 
in the SK detector in the following. 

Xf = x Q , (D-ll) 
Vf = 2/0 + £T sk sin v (p) (D-12) 
Zf = z Q + £T sk cos0„(i>). (DT3) 

If we carry out the Monte Carlo Simulation, following 
the flow chart in Fig. 1171 then, we obtain N event , the num- 
ber of the neutrino events generasted in the SK detector. 
The ratio of the selected events to the total trial is given 
as 

-Rmontc (COS = A^ even t (COS Q v (v) ) /-^sample (COS 

(D-14) 

Comparison between i?thcor and i? m0 nte hi Table 3 shows 
that our Monte Carlo procedure is valid. 



Entry 




N^NsamplB, 

YES 
Determination of 
Point in the plane ABCD 
by using I 1 and I 2 




Determination of Tweigtu 
p iiv. f p and £ sk 




Fig. 17. Flow Chart for the determination of the interaction 
points of the neutrino events inside the detector 
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